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Abstract A joint Mexican, Russian, and American research effort has been initiated to develop 
new methods to remotely sense atmospheric parameters using ground-based, aircraft, and satel- 
lite observations. As a first step in this program, ground-based spectrophotometric measurements 
of the direct solar radiation have been obtained for the extremely polluted Mexico City atmos- 
phere for the period of April-June 1992. These observations were made at more than 1300 chan- 
nels in the spectral range of 0.35-0.95 pm. In the ultraviolet (UV) portions of the spectrum (e.g., 
0.35 pm), aerosol optical thicknesses were found to range between 0.6 and 1.2; in the visible 
portion of the spectrum ( e g., 0.5 pm) they ranged from 0.5 to 0.8; and in the near-infrared 
(NIR) spectra (e.g., 0.85 pm), values of 0.3 - 0.5 were found. Applying a spectral optical depth 
(SOD) model of x(^) = C + A values of 1.55 < oc < 1.85 were obtained for polluted, cloud- 
less days, with values of 1.25 < <x < 1.60 on days with haze. The aerosol particles in the polluted 
Mexico City atmosphere were found to be strongly absorbing, with a single-scattering albedo of 
0.7 - 0.9 in the UV, 0.6 - 0.8 in the visible portion of the spectrum, and 0.4 - 0.7 in the NIR 
These values are possibly consistent with a high soot concentration, contributed both by vehicular 
traffic and heavy industiy. Analysis of the measured aerosol SOD using the optical parameters of 
an urban aerosol model permits the concentration of aerosol particles to be estimated in the verti- 
cal column; a maximum value of 3 x 10 ^ cm"^ was found. This concentration of aerosol parti- 
cles exceeds that found in most other regions of the globe by at least an order of magnitude. Near 
the ground the aerosol size distributions measured using an optical particle counter were found to 
be strongly multimodal. 


1. Introduction 

Aerosols, both anthropogenic and natural, are considered to 
be uncertain but potentially important climate feedback 
mechanisms [Rasool, 1973; Schneider , 1976; Hansen and 
Lacis , 1990; Mitchell, 1991]. Aerosols reflect sunlight back 
toward space, thereby cooling the Earth's surface. They also 
cool the planet indirectly by providing cloud condensation 
nuclei (CCN) upon which water vapor condenses into cloud 
droplets. Large concentrations of small droplets make the 
clouds more reflective to sunlight [Coakley et al , 1983]. 
Furthermore, observations demonstrate that aerosols increase 
cloud reflectivity by inhibiting rainfall and prolonging cloud 
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lifetimes [Albrecht, 1989]. On the other hand, some aerosols 
may strongly absorb solar radiation, thereby heating the terres- 
trial atmosphere [Kondratyev et al, 1974a, b, 1977, 1981; 
Welch et al., 1981; Vasilyev et al, 1993]. Such aerosol parti- 
cles located on the surface of water droplets also may increase 
the absorbtivity of cloud layers [Twohy et al, 1989], The net 
result is that changing amounts of anthropogenic aerosols may 
contribute to changed radiative forcing and thus to climate 
change. 

Current climate forcing by anthropogenic sulfate is esti- 
mated to be of the order of -1 to -2 W/m^, which is compara- 
ble in magnitude but opposite in sign to greenhouse gas forc- 
ing [Charlson et al, 1992]. Aerosol forcing may make a 
significant contribution to offsetting global greenhouse 
wanning; however, there are large uncertainties in our under- 
standing of the magnitude of these effects. 

The radiative balance in the atmosphere is strongly influ- 
enced by the scattering and absorbing characteristics of aero- 
sols. Nonabsorbing aerosols only increase the albedo of the 
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surface-atmosphere system, thereby cooling the surface 
[Kondratyev et al, 1973]. However, for absorbing aerosols 
and especially highly absorbing aerosols such as soot the ef- 
fects become more complicated. Indeed, the radiative balance 
of the surface-atmosphere system depends upon both the sign 
and the magnitude of the difference between the surface al- 
bedo and the reflectivity of the aerosol layer [Kondratyev et 
al, 1973, 1977; Kondratyev and Vasilyev, 1975]. For cases of 
a dark underlying surface (i.e., low albedo), the presence of 
aerosols produces a net increase in the atmosphere-surface 
system albedo. On the other hand, for bright surfaces, some 
aerosols may produce a net decrease in system albedo. 

A long series of airborne observations of radiative flux di- 
vergences have been carried out in Russia for a large variety 
of atmospheric conditions [Kondratyev et al. y 1974a, b, 1977, 
1981; Vasilyev , 1993]. For very dusty or highly polluted at- 
mospheres the solar radiation absorbed by aerosols may reach 
values equivalent to gaseous absorption. Similarly, aerosol 
scattering of solar radiation may reach values equivalent to 
scattering by the molecular atmosphere. However, the meas- 
urements show large geographical and seasonal variations of 
aerosol effects and hence of aerosol radiative forcings on cli- 
mate. 

Atmospheric aerosol particles generated from dust storms 
[ Carlson , 1979] and smoke from forest fires [Chung and Le , 
1984] are important regional climatic variables. Large parti- 
cle concentrations produced from these types of events often 
are found with optical thicknesses greater than one. These 
aerosol particles also are transported across great distances 
from the source [Fraser et al., 1984], and they perturb the ra- 
diative balance significantly, both on local and on regional 
scales. Heavy dust loading can change the top-of-the- 
atmosphere clear sky shortwave and longwave radiative exis- 
tence by 40-90 W/m r and 5-20 W/m 2 , respectively [Ackerman 
and Chung , 1992]. 

Sulfate is one of the important sources of tropospheric aero- 
sols. Sulfur dioxide emissions are produced by both anthro- 
pogenic [Schwartz, 1988] and by natural [Chari son et al . , 
1987] sources. Schwartz [1988] argues that anthropogenic 
SC >2 emissions cannot be neglected in analyzing the global 
cloud albedo-climate issue [Charlson etal. , 1992]. 

Large urban areas often are sources of extremely heavy 
pollution, affecting both the radiation balance and the climate 
for extended regions around them and the health of the popu- 
lation within them [Ensor et al . , 1974]. For example, the 
aerosol concentrations in the Grand Canyon are contributed to 
by both the Los Angeles basin, hundreds of miles to the west, 
and by Phoenix to the south. Pollution from North America is 
routinely tracked all the way to Europe, and both North 
America and Europe strongly contribute to Arctic haze. Urban 
pollution sources are continuous, thereby contributing to the 
overall regional background aerosol concentrations. 

The Mexico City atmosphere is widely held to be the most 
heavily polluted region on the planet. The source of the aero- 
sols is a combination of heavy vehicular traffic and heavy in- 
dustry, neither of which has any significant emission require- 
ments. Indeed, the heavy pollution found in Mexican cities is 
an international issue and a consideration in the North Ameri- 
can Free Trade Agreement (NAFTA). Therefore the first site 
for a joint Mexican, Russian, and American research effort, 
initiated to develop new methods to remotely sense atmos- 
pheric parameters using ground-based, aircraft, and satellite 


observations, is Mexico City. As a first step in this research 
program, ground-based spectrophotometric measurements of 
the direct solar radiation have been obtained in Mexico City 
for the period of April- June 1992. 

Note that most retrieved measurements of aerosol concen- 
trations and size distributions reported in the literature are 
based on single or, at most, a few wavelengths. In contrast, 
the present observations have been made at more than 1300 
channels in the spectral range of 0.35 - 0.95 pm. This wide 
range of wavelengths allows a significant improvement in 
characterizing the polluted atmosphere. Section 2 discusses 
the instrumentation, and section 3 describes the sites and 
conditions of the observations. Section 4 outlines the meth- 
odology; section 5 reports the results; and conclusions are pre- 
sented in section 6. 

2. Instrumentation 

Measurements of the spectral optical depth (SOD) in Mex- 
ico City during April- June 1992 were made using the spectral 
information measuring system (SIMS), consisting of (1) an 
optical unit, including a projection optical system and an obtu- 
rator of the light beam; (2) a grating monochromator, (3) an 
electronic amplifying and directing unit; and (4) a registering 
and preliminary processing unit (i.e., PC computer). The 
principal element of SIMS is the grating monochromator (the 
Russian K-2 or K-3), which was designed using the vertically 
symmetric scheme of Ebert-Fastie. The main parameters of 
the K-2 and K-3 instruments are listed in Table 1 . Schematic 
and interior views of the K-3 monochromator are shown in 
Figure 1. The K-2 and K-3 spectrophotometers were designed 
and manufactured in the Laboratory of Experimental Atmos- 
pheric Optics at the Institute of Physics at St. Petersburg State 
University (IPh SPbSU) [Vasilyev etal, 1987]. 

The two concave spherical mirrors have focal lengths of 193 
inches (482.5 mm) and the plane grating has 600 lines/mm. 
The first order of diffraction is used. With a spectral resolu- 
tion of 0.005 (i.e., 1/200) the output slit linear resolution is 33 
A/mm. By rotating the grating, output from the various 
wavelengths is measured. The measuring time is 7 s, and the 
time of returning the grating to its starting position is 3 s. 
Radiation dispersed by the grating falls on three output slits 
and is received by three photomultipliers (the UV, visible 
(VIS), and NIR portions of the spectrum). Each output slit 
has a glass correction filter to "cut off" the higher orders of 
diffraction. 

Measurements of radiative fluxes are made using opaque 
glass mounted on the entrance objective. Their "cosine law" 
was found to deviate by less than 0.5% for angles less than 
70°. Measurements of direct solar radiation then were made 
using this opaque glass and a long tube with an entrance dia- 
phragm to limit the field of view of the instrument. 

During the measurement process, instrument stability was 
controlled using a lamp mounted inside the monochromator. 
Before the observations were taken, this calibration lamp was 
compared to a standard source Russian SI-8 lamp. In turn, 
this standard lamp was calibrated using a "blackbody" source 
of light at the Shtemberg Astronomical Observatory at Mos- 
cow State University. 

The "warm-up" time of the instrument is about 1/2 hour. 
During the measurements, variations of the signal were found 
to be 1-2%. The lamp inside the monochromator was not 
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Table 1 Main Parameters of K-2 and K-3 Instruments 


Instrument Parameter 

Value 

Dispersing element 

Diffraction Grating, 600 mm’ 1 

Order of diffraction 

1 

Spectral range 

0.35 - 0.95 pm 

Sensors 


three photomultipliers 

UV, 0.35 - 0.42 pm 

for three channels 

VIS, 0.42 - 0.70 pm 


IR, 0.70 - 0.95 pm 

Spectral resolution 

1 : 200 

Discrimination Sensitivity 

0. 1 mW/cm 2 pm 

Field of view 

1 0 x 10’ or 180° 

Measuring time 

7s 

Operating voltage 

27 V 

Consumed power 

200 W 


thermostated; since luminosity is proportional to T 4 , this 1- 
2% signal variability may be traced to a 1/4- 1/2% variability 
of the temperature (T) of the lamp’s filament. In absolute val- 
ues this variation represents changes of the lamp's filament 
temperature of 6°-12°C, referenced to the original temperature 
of about 2400°C. As a result, all measurements taken using 
SIMS had an accuracy of about 1-2%. The calculated values 
of SOD have a similar accuracy. 

3, Measurement Sites and Conditions 
of Observations 

The north side of Mexico City is a region of heavy industry 
and the main source of industrial pollution. The main obser- 
vational site for the measurements is the Observatory of Solar 
Radiation (ORS) at the Institute of Geophysics (IGf) of the 
Mexican Autonomous National University (UNAM), located 
in the southern part of Mexico City. This so-called 
"ecological zone" is a region in which air pollution is lower 
than in the center of the city. Observations were conducted 
from the roof of the Institute of Geophysics, UNAM, and in 
the central part of Mexico City at Tezozomoc Park, in one of 
the most polluted regions of the city. 

Another series of measurements was taken at Tlamacas Sta- 
tion (4 km above sea level) on the slope at the Popocatepetl 
volcano. These observations were designed both for retriev- 
ing the background optical thickness measurements and for 
obtaining an absolute calibration of the K-3 spectrophotome- 
ter, which is necessary if the observations are conducted using 
the so-called “short method” (see section 4). However, the 
choice of the Tlamacas Station site proved to be unfortunate 
because measurements of the Sun had to be made in the di- 
rection of the city of Pueblo which is the second largest indus- 
trial center of Mexico. The air above this region also was 
found to be highly polluted and unusable for background re- 
trievals and calibration purposes. 

A series of observations also were taken in the presence of 
high cirrus cloudiness; therefore some estimates of cirrus SOD 
are included in the results. 


The meteorological conditions prevailing during the period 
April-June 1992 generally permitted only a morning series of 
measurements, with the one exception being the afternoon of 
June 14. During a 2-month observational period a series of 24 
observations were made, or approximately one series of ob- 
servations each 3 days. On the other hand, there were periods 
of cloudiness lasting for up to 10 days when no measurements 
were made. 

4. Methodology 

The "long method" for the determination of spectral optical 
depth was used in this investigation [e.g,, Liou, 1980]; that is, 
the direct solar radiation was measured at a series of different 
solar zenith angles, corresponding to different "air masses" 
(i.e., optical path lengths). Using Beer's law, the observed 
fluxes of direct solar radiation are expressed as 

S(X) = S 0 (X) e- 1 ^) m ( z ) > ( 1 ) 

where X is wavelength; S 0 (X) is the extra-atmospheric solar 
spectral flux (i.e., the so-called "spectral solar constant"), t(X) 
is the spectral optical depth, z is the solar zenith angle, and 
m(z) is the air mass. Taking the logarithm of this equation. 

In S(X) = In S Q (X) - t(X) m(z) , (2) 

a linear regression is obtained which may be solved using the 
least squares method. 

The advantage of the “long method” is that the spectral op- 
tical depths (SODs) are obtained without requiring any cali- 
bration of the instrument, because they are relative, not abso- 
lute measurements. The linearity of the instrument was 
carefully examined in laboratoiy studies, and the stability of 
the instrument was controlled during the observations using a 
lamp mounted inside the monochromator (see section 2). The 
calibration of the instrument is critical for observations using 
the “short method,” in which values of SOD are 
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Figure 1. (a-c) Schematic and (d) interior view of the K-3 monochromator, (b) Numbers 1-3, photomul- 
tipliers; 4, quartz windows; 5, correction filters, (c) Number 1, opaque glass; 2, obturator; 3, income slit; 
4, solar blend; 5 and 9, inner diaphragms; 6 and 7, opaque screens, 8, objective; 20, standard lamp. 
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S(X) 


(3) 


Figure 2 shows an example of the Langley plots in the ultra- 
violet (UV), visible (VIS), and near-infrared (NIR) portion of 
the spectra measured in Mexico City on April 12, 1993. 

Most of the observations were made in the morning from 
0700 to 0900 Local Standard Time (LST), during which time 


the atmosphere generally is expected to be stable. Neverthe- 
less, the optical state of the atmosphere does begin to change 
during this period of time as the industrial plants begin op- 
eration and as the traffic density increases. Such variations 
may lead to errors in the retrieved results using the "long 
method." The optical path length varied from about 5 to about 
2, during which time the solar elevation ranged from 5° to 
30°. Because of the low latitude of Mexico City the Sun rises 
at a large angle with respect to the horizon; hence the change 
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of the Sun's azimuth does not exceed 20° during this meas- 
urement time. As a result, there is no significant azimuthal 
influence on the observational data. 

The spectral optical depth x(X) is composed of four parts: 
(1) molecular (Rayleigh) scattering, (2) molecular absorption, 
(3) aerosol scattering, and (4) aerosol absorption. If one sub- 
tracts the molecular (Rayleigh) scattering part from t(l), the 
remainder, the so-called “reduced SOD,” characterizes aero- 
sol attenuation (scattering and absorption) and molecular ab- 
sorption. 

It is common to apply the Angstrom turbidity formula 
[Angstrdm, 1961; Iqbal , 1983] to the analysis of aerosol 
measurements: 


^) = PX* QC , (4) 

where P is the so-called Angstrom turbidity coefficient (note 
that this coefficient is equal to t(X= 1 pm)). Angstrdm [ 1961 ] 
suggested oc = 1.3±0.5asan average value and showed that 
if a Junge size distribution N(r) of aerosol particles is as- 
sumed, 


N(r) = N 0 r‘ v , (5) 

where r is the particle radii and N 0 and v are constants, then 
an exact relationship between x in equation (4) and v in 
equation (5) exists, namely, 

v = a + 2 . (6)‘ 

Therefore assuming that the aerosol size distribution is 
Jungian, the slope of the averaged size distribution can be es- 
timated from the spectral dependence of the optical thickness 
measurements. 

Unfortunately, this Angstrflm formula approach does not 
represent the observations well, because it is based only on 
Mie scattering by the large particles in the size distribution. 
Other problems which limit the applicability of this approach 
are: (1) the absorption of light by aerosol particles is not ac- 
counted for, (2) most size distributions are multimodal; and 
(3) giant particles (i.e., strong haze, fog) tend to scatter light 
almost neutrally (i.e., wavelength independent). As a result, 
an alternative formula has been suggested [e.g., Shifrin and 
Minin , 1957] of the form 


tn 




t(X) = C + AX'* , (7) 

where C and A are constants. Using this approach, the Ang- 
strom turbidity coefficient, computed as b=t(l) at X = 1 pm, is 


10 4 r , — r 

1.50 2.00 2.50 3.00 3.50 4.00 

m(Z ) 

Figure 2 . Examples of “Langley plots” (or “Bouguer lines”) 
taken m Mexico City on April 12, 1992, for (a) UV at 1 = 377 
nm; (b) visible at 1 = 545 rnn; and (c) NIR at 1 - 814 nm 


P = C + A (8) 

In equation (7), C is the contribution due either to absorp- 
tion by all aerosol particles or to scattering by giant particles. 
Assuming C to be due only to absorption, the spectral single- 
scattering albedo TU 0 (AJ may be estimated for aerosol par- 
ticles as 

c 

t(A) 



(9) 
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Figure 3. Aerosol Spectral Optical Depth (SOD) in Mexico City in April 1992. Number 1, April 14; 2, 
Ami 12 3 April 16 4, April 15, 5, WMO Standard Radiation Atmosphere (1986), 6, SOD in the Elter- 
man model (1968) With the visibility of S m = 25 km; SOD measured at different sites in Russia 
[Kondratyev et al., 1965; Kondratyev ed„ 1969] at different heights above sea level, 7, 45 m; 9, 3410 m, 
10, 3700 m, and 8, SOD measured in Tucson, Arizona [Box et al., 1981] in May 1977. 


If scattering by very large aerosol particles cannot be ne- 
glected, then the estimate in equation (9) will be reduced; in 
such cases, equation (9) is considered as a minimum estimate. 

5. Results 

5.1. Aerosol Spectral Optical Depths 

Reduced spectral optical depth (total spectral optical depths 
minus Rayleigh depths) observations taken in Mexico City on 
April 12, 14, 15, and 16 are shown in Figure 3 (curves 14). 
Absorption bands of ozone, water vapor, and oxygen are rec- 
ognized. The spectral signature in the UV-blue spectral range 
probably is caused by the absorption by nitrogen dioxide, of- 
ten present in large concentrations in urban regions. Values of 
aerosol SOD measured in Mexico City are much larger than 


those found in most other regions of the world [e g., Kondrat- 
yev et al , 1965, 1974a, b; Box et al., 1981]. Curves labeled 
5-10 in Figure 3 are used as reference spectra of typical back- 
ground values found in other regions of the globe. Curves 7, 
9, and 10 were obtained in unpolluted regions of Russia at 45 
m, 3410 m, and 3700 m above sea level. Curve 8 was meas- 
ured in Tucson, Arizona, in May 1977 [Box et al 1981]. Ob- 
servations made in different years in Tucson have similar val- 
ues [Shaw etal, 1973; King etal , 1978; King, 1982]. Curve 
5 shows continental aerosol SOD values according to the 
Standard Radiation Atmosphere [World Meteorological Or- 
ganization (WMO), 1986]; and, finally, curve 6 represents the 
results of SOD calculations made using the polluted aerosol 
model of Elterman [1968] which is similar to the Q.S. Stan- 
dard Atmosphere. This model assumes that the aerosol con- 
centration corresponds to a visibility of 25 km, which is con- 
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Table 2 Solution of equation (7) Using the Method of Variation 
of Parameter a 2 


Olz 

c 2 

A2 

r 2 

S 2 

0.10 

-4.524 

4.790 

0.9728 

0.0948 

0.20 

-1.987 

2.259 

0.9758 

0.0845 

0.30 

-1.141 

1.419 

0.9789 

0.0739 

0.40 

-0.719 

1.003 

0.9819 

0.0636 

0.50 

-0.466 

0.755 

0.9845 

0.0546 

0.60 

- 0.297 

0.592 

0.9869 

0.0464 

0.70 

-0.176 

0.477 

0.9891 

0.0386 

0.80 

-0.086 

0.392 

0.9910 

0.0319 

0.90 

-0.016 

0.327 

0.9926 

0.0263 

1.00 

0.041 

0.276 

0.9940 

0.0216 

1.10 

0.087 

0.235 

0.9950 

0.0180 

1.20 

0.125 

0.201 

0.9958 

0.0151 

1.30 

0.158 

0.174 

0.9962 

0.0136 

1.40 

0.186 

0.151 

0.9965 

0.0127 

1.50 

0.210 

0.132 

0.9963 

0.0132 

1.60 

0.232 

0.115 

0.9960 

0.0145 

1.70 

0.251 

0.101 

0.9952 

0.0172 

1.80 

0.268 

0.089 

0.9943 

0.0206 

1.90 

0.283 

0.079 

0.9930 

0.0251 

2.00 

0.297 

0.069 

0.9915 

0.0306 

2.10 

0.309 

0.061 

0.9898 

0.0369 

2.20 

0.321 

0.055 

0.9877 

0.0443 

2.30 

0.331 

0.048 

0.9854 

0.0525 

2.40 

0.341 

0.043 

0.9829 

0.0616 

2.50 

0.350 

0.038 

0.9801 

0.0715 

2.60 

0.358 

0.034 

0.9771 

0.0823 

2.70 

0.366 

0.031 

0.9738 

0.0937 

2.80 

0.374 

0.027 

0.9704 

0.1059 

2.90 

0.381 

0.024 

0.9668 

0.1186 

3.00 

0.387 

0.022 

0.9629 

0.1321 

3.10 

0.394 

0.019 

0.9589 

0.1461 

3.20 

0.399 

0.017 

0.9547 

0.1606 

3.30 

0.405 

0.016 

0.9504 

0.1757 

3.40 

0.410 

0.014 

0.9459 

0.1911 

3.50 

0.416 

0.013 

0.9412 

0.2070 

3.60 

0.420 

0.011 

0.9365 

0.2231 

3.70 

0.425 

0.010 

0.9316 

0.2397 

3.80 

0.430 

0.009 

0.9265 

0.2566 

3.90 

0.434 

0.008 

0.9214 

0.2736 

4.00 

0.438 

0.007 

0.9163 

0.2908 


sidercd to be moderate pollution. The values of SOD of the 
polluted atmosphere in Mexico City are significantly larger, 
representing visibilities of about 3 km. 

Note the spectral dependence of the curves shown in Figure 
3. As the optical thickness increases, the slopes of the curves 
likewise increase. Aerosol absorption typically is large for the 
type of particles typically found in the polluted urban envi- 


ronment. Therefore when the aerosol concentrations are low, 
the role of aerosol absorption in the attenuation of light domi- 
nates over scattering by aerosol ensembles. Then the spectral 
dependence of aerosol optical depth follows closely the de- 
pendence of the imaginary part of the aerosol refractive index. 
For many aerosols (salts, soils, and soot) this spectral depend- 
ence is nearly neutral in the 0.4-to 1.0-pm spectral range. In 
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Figure 4. Comparison of different approximations of SOD (explanations are given in the text). 


such cases, the absolute values of SOD generally do not ex- were out of the confidence level of 0.995 also were rejected, 

ceed t = 0. 1-0.2. Such is the case for curves 6-9 in Figure 3. The Fisher F-criteria were applied to the ratio of interseries 

For curves 5 and 6, with an optical thickness of about 0.2- deviation to a random one as a fmal rejection test (confidence 

0.55, scattering becomes important; hence the spectral de- level of 0.995). As a result, half of the measurements were 
pendence becomes nonneutral. As seen in curves 1-6, the rejected. This should be expected, because the optical prop- 

spectral signature is due to the increasing role of scattering of erties of the morning Mexico City atmosphere often do not 

light in optically thicker aerosol layers, especially for wave- remain stable for extended periods of time, 
lengths in the range of 0.35-0.70 pm. In the highly polluted 

Mexico City atmosphere the amount of radiation attenuated by 5 , 2 , Spectral Dependencies of Measured Optical Depth 
aerosols is of the same order of magnitude as that scattered by 

the molecular atmosphere (i.e., Rayleigh scattering). The attenuation of light by aerosol particles is caused by 

As noted above, during the processing of the measurements, both scattering and absorption processes. Usually it is im- 

so-called “Langley plots" (or "Bouguer lines") were computed possible to separate one effect from the other, hence Angstrdm 

(see Figure 2). These are regression lines of the logarithm of [1961] suggested the use of the Angstrdm turbidity formula, 

measured intensity against optical path length (i.e., air mass). given by equation (4). This equation expresses the spectral 

The correlation coefficient obtained from this analysis was dependence of optical depth versus wavelength only in cases 

found to be greater than 0.95, indicating a relatively low level of weak air pollution (refer to the parameters of this equation 

of experimental error. However, one obtains reliable values as aj and Cj = P). For the strongly polluted atmosphere of 

of SOD only when the optical parameters of the Earth's at- Mexico City it is necessary to use equation (7) (refer to these 

mosphere do not change during the observational period. three parameters as a 2 , C 2 , and A 2 ). Note that this equation 

Therefore the measurements were examined by applying sta- is nonlinear. To determine all three coefficients simultane- 

tistical processing methods; for example, data were rejected ously, one of these parameters is fixed and the other two are 

with large standard deviations. Using the Pierson x 2 -criteria computed. It is usually most convenient to fix the power a in 

(confidence level of 0.995), all measurements with large er- equation (7), since it is known that this parameter may vary 

rors in determination of the SOD were rejected. Rejected only between 0 and 4. The calculations are performed in the 

measurements were those which had only a few observations following way: Using the least squares method, the solution 

or those in which there were only a few measurements at dif- is found for fixed values of a, varying from 0 to 4, with the 

ferent solar zenith angles (i.e., limited change of optical path step size equal to 0.1. That solution is chosen which has the 

length intervals). Then, using the Student t-criteria, meas- highest correlation coefficient (or the lowest standard devia- 

urements whose extrapolated direct solar radiation fluxes tion). Of course, differences between the alternative solutions 
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Table 3 Results of the Approximation of SOD Taken in Mexico City 
in 1992 


Angstrftm Formula Formula (7), Variation of a Formula (b), Determination of 


Date, 

a i 

Cl 

s. 

<*2 

c 2 

^2 

s 2 

03 

C3 

A 3 

^3 

Morning 


xlO 

xlO 2 


xlO 

xlO 

xlO 2 


xlO 

xlO 

xlO 2 

Series 












May 01 

1.25 

3.12 

9.5 

1.84 

1.96 

1.45 

4.8 

1.76 

1.83 

1.59 

5.0 

May 16 

1.15 

2.95 

4.1 

1.54 

1.62 

1.64 

2.8 

1.69 

1.85 

1.39 

3.0 

May 21 

1.05 

2.85 

4.7 

1.73 

2.14 

1.07 

2.3 

1.57 

1.95 

1.24 

2.5 

May 28 

1.40 

1.94 

5.3 

1.78 

1.06 

1.16 

3.9 

1.86 

1.12 

1.06 

4.0 

Average 

1.22 

2.69 

4.7 

1.75 

1.71 

1.29 

1.8 

1.72 

1.69 

1.32 

3.8 

(a) 












April 12 

0.87 

3.30 

3.7 

1.06 

1.06 

2.35 

3.4 

1.57 

2.03 

1.19 

4.5 

April 14 

0.88 

4.00 

5.6 

0.88 

0 

4.00 

5.6 

1.57 

2.70 

1.49 

12.3 

May 25 

1.02 

2.30 

4.2 

1.23 

0.77 

1.66 

4.0 

1.61 

1.66 

0.96 

4.6 

Average 

0.90 

3.23 

11.0 

0.90 

0 

3.23 

1.1 

1.59 

2.13 

1.21 

7.1 

(b) 












May 15 

0.69 

3.40 

14.4 


solution 

not found 


1.28 

2.91 

1.07 

12.0 

May 27 

0.67 

3.60 

8.9 

0.67 

0 

3.60 

8.9 

1.61 

2.22 

1.75 

18.0 

June 14 

0.60 

3.23 

3.5 


solution 

not found 


1.34 

1.52 

0.78 

1.5 

Average 

0.60 

3.57 

2.4 

0.80 

1.35 

2.31 

2.3 

1.41 

2.22 

1.20 

10.5 

(c) 












Average 

0.92 

3.12 

2.5 

1.38 

1.81 

1.55 

1.3 

1.45 

1.98 

1.25 

1.3 

(a-c) 












April 25 

0.16 

18.5 

17.3 


solution 

not found 


0.72 

19.6 

1.43 

14.8 

May 02 

0.04 

8.45 

23.7 


solution not found 


1.05 

4.81 

0.80 

28.6 

May 13 

0.13 

9.69 

13.8 


solution 

not found 


0.65 

6.20 

1.72 

13.8 

Average 

0.13 

12.2 

_ 


solution 

not found 


0.81 

10.2 

1.32 

19.1 

m 












st be statistically significant. Only 

in cases 

of strong pol- 

mental data, a value of a = 

0.925 would have resulted 


lution is it statistically reliable to apply this approach. 

Table 2 illustrates how the method is applied in practice. 
The data used are the average of observations taken in Mexico 
City during the whole observational period of 1992. In this 
case, the maximum correlation coefficient R2 corresponds to 
the value of a = 1 .4 with the minimum standard deviation S2 = 
0.0127, To achieve higher accuracy, the calculations may be 
repeated in a narrower range of a with the step size equal to 
0.01. Note that a solution with C2 = 0 corresponds to the 
Angstrom turbidity formula, equation (4). Such a solution 
would correspond to the value of a = 0.925, with the value of 
Angstrom turbidity coefficient A2 - 0.31 and the standard 
deviation of 0.025, which is twice as large as in the previous 
case. This means that had the Angstrom turbidity formula 
(equation (4)) been used for the approximation of the experi- 


of the correct value of a = 1.4. Note also that the Angstrom 
turbidity coefficient is equal to 0.31; this value is about 10% 
low, because equation (7) yields C2 + A2 = 0.34. However, 
once again note that the approach of varying parameter 012 is 
applicable only in cases of strong air pollution (when t(^=0.5 
pm) >0.6). 

For cases of weaker air pollution, another approach may be 
used. One may assume an exponential form of t(l): 

a fX 

x(X) = c 3 e 3 (io) 

This equation is linearized taking its logarithm: 

In t(l) = In C3 + a3/l (11) 

In this case, the parameter C3 and a3 are determined using the 
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Table 4: Characteristics of the Atmospheric Turbidity Q Mexico City in April- June 
1992 


Groups of Days 

Cj (Angstrfim 
Formula) 

A r + C* 

A, + 

0.5xt(A= 0.5 pm) 

Cloudless days, 1 gr. (a) 

0.27 

0.30 

0.30 

0.31 

Cloudless days, 2 gr. (b) 

0.32 

0.32 

0.33 

0.31 

Days with haze (c) 

0.36 

0.37 

0.34 

0.27 

All cloudless days (a-c) 

0.31 

0.34 

0.32 

0.30 


least squares method. Then if the value of C 3 is known, it is 
easy to linearize equation (7) by taking its logarithm; in par- 
ticular, 

ln[T(X)-C 3 ] = lnA 3 -a 3 lnX (12) 

Parameters A 3 and 013 are computed using the leastsquares 
method. Mathematically, this means that the approximation 
is based upon equation (7), but it is restricted to pass through 
the point (0, C 3 ) on the plot of In i(k) versus 1 /X (see equa- 
tion ( 11 )). 


Figure 4 illustrates the errors which occur using this ap- 
proach. Again we use the averaged values of x(X) calculated 
for all of the series of measurements taken in Mexico City in 
1992. This averaged curve and its extension to 1/k -> 0 , us- 
ing the method of variation of a , is shown by squares in Fig- 
ure 4. The exponential approximation for the spectral range 
of observations (0.35 - 0.95 Jim) and its corresponding exten- 
sion to 1/A. — ► 0 is shown by the solid line. The approxima- 
tion using equation (7) with the constant C 3 , determined us- 
ing equation (1 1), is shown by the crosses in Figure 4. Within 
the spectral range of the measurements, differences between 



C 


Figure 5. Dependence of the parameter a versus the parameter C: plus, clear days, case “a”; circle, clear 
days, case “b”; square, days with haze; and triangle, days with cirrus clouds. 
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Figure 6. Aerosol SOD observed in Mexico City in 1992 
divided into groups according to the values of the parameter a 
(the order of the dates from the top to the bottom: (a) May 1, 
16, 21, 28, and Rayleigh; (b) April 14, 12, and May 25; (c) 
May 13, 2, 27, 15, and June 14 (morning series). 


these approaches are not statistically significant; the curves ef- 
fectively coincide. The difference between C 2 = 0.181 and 
C 3 = 0.198 is about 9%, and the value of a varies from 1.40 
to 1.45, about a 3% difference. Taking into account the fact 
that the observational errors are of the order of several per- 
cent, these differences may be considered as insignificant. 

The results of the different approximations of the aerosol 
spectral optical depths measured in Mexico City are given in 
Table 3. In all cases the Angstrom turbidity formula provides 
low values of both the parameter a and the Angstrom turbid- 
ity coefficient. The method of variation of the parameter a 
provides a solution only in cases of high air pollution, while 


the approach of determining the parameter C 3 from the expo- 
nential approximation works in all cases under consideration. 

For various degrees of atmospheric turbidity. Angstrom de- 
fined the following values of turbidity coefficient: clean, « 
0.0; clear, » 0.1; turbid, « 0.2; and very turbid, * 0.4. The 
corresponding visibilities are 340 km, 28 km, 11 km, and <5 
km, respectively. Table 4 gives the value of the turbidity co- 
efficient determined using the different approaches. Values of 
turbidity coefficient range from 0.27 to 0.37, with visibilities 
in the range of about 7 km to 13 km, corresponding to a very 
turbid atmosphere. Visibilities were observed to be about 3 
km or less, which coincide (within the limits of experimental 
error) with the estimate made using Angstrom turbidity coef- 
ficient. 

The parameters CI 3 and C 3 are related, as shown in Figure 
5. The solid line is a power law which is computed using the 
least squares method. 

The spectrophotometric measurements may be divided into 
four groups. Some polluted, cloudless days are found to have 
a strong dependence of x versus X, with values ranging from 
1.70 < a < 1.85. Other polluted, cloudless days are found 
with a weaker dependence of t versus X, with 1.55 < a < 
1.65. For days with haze, there is a still weaker dependence 
of t versus X, with values of 1.25 < a < 1.60. Finally, when 
cirrus clouds are present, there is nearly neutral dependence 
of x versus X, with values of 0.65 < a < 1.05. These differ- 
ences in the values of a may be caused by several reasons. 
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Figure 7. Aerosol SOD observed in Mexico City in 1992: (a) 
averaged according groupings in Figure 6; (b) averaged over 
all the observations on clear days. 
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Figure 8. SOD of cirrus clouds observed in Mexico City on April 24, 1992: number 1, observed SOD, 2, 
average curve for all cloudless days (Figure 7d); and 3 SOD of only cirrus clouds. 
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Figure 9. Temporal variations of aerosol SOD at three wavelengths: 350 nm, 500 nm, and 850 nm for the 
period beginning April 10, 1992. 
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Figure 10. Comparison of SOD measured in Mexico City on June 14 in the morning (top curve) and in 
the evening (bottom curve). 


First, for high relative humidity, the aerosol particles may be- 
come covered by water, decreasing the absorption but increas- 
ing scattering in the forward direction (e g., detected as the 
presence of a haze). Second, the aerosol particles may coagu- 
late, resulting in a variation in size distribution. Increasing 



the numbers of large particles will produce a decrease in the 
slope of the x versus k dependence. Another factor also needs 
to be considered. Multiple scattering increases with increas- 
ing particle concentration (i.e., increasing t), changing the 
slope of the x versus X dependence. All of these factors may 



Figure 11. (a) Comparison of aerosol SOD measured on May 28: number 1 in the center of Mexico City; 
and 2 at the territory of Cd. Universitaria; (b) the ratio of intensities of direct solar radiation. 
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Figure 12. Calculated aerosol SOD for different amounts of aerosol particles (0.5, 1.8, 3.0, and 5.0 x 10' 9 
cm* 2 ) in the vertical column of the atmosphere with a crosssection equal to 1 cm . 


contribute in the Mexico City atmosphere. The process of co- 
agulation leads to a decrease in particle concentration, a slight 
improvement in atmospheric transparency and a decrease in 
the slope of the x versus X dependence. On the other hand, an 
increase in particle concentration leads to a deterioration of 
atmospheric transparency and an increase in the slope of the x 
versus X dependence, due to multiple scattering. 

Figure 6 shows the results of the measurements of SOD di- 
vided into each of four groups of days mentioned above. As a 
point of reference, the bottom curve in Figure 6a represents 
the Rayleigh scattering contribution for the Mexico City at- 
mosphere; note that the effects of pollution upon optical 
thickness far outweigh those from Rayleigh scattering for all 
observations and at all wavelengths. In Figures 6a and 6b, 
note that the optical thickness is approximately equal at the 
larger wavelengths but varies substantially at the shorter 
wavelengths. According to Angstrom's law this means that 
stronger episodes of pollution are characterized by steeper 
slopes of the aerosol particle size distribution. As noted 
above, in Figure 6a, either there were larger concentrations of 
relatively small particles in the atmosphere than in Figure 6b 
or the slope of the dependence x versus X changed due to the 
increased role of multiple scattering; or probably, both effects 
were present. The conclusion is that strong pollution is char- 
acterized by increasing concentrations of small particles and 
by the increasing role of multiple scattering. As pollution 
weakens (e g.. Figure 6c), the hypothesis is that the concen- 
tration of smaller-size particles decreases with respect to the 
larger-sized particle. In this case, the spectral dependence of 
SOD flattens out, being influenced by the giant particles. 
When cirrus clouds are present, the spectral dependence of 
SOD becomes almost neutral. 

Figure 7 shows averaged curves for the groups shown in 
Figures 6a-6c and an average curve for all the cloudless days. 
Figure 7a shows that the air pollution decreases from case a to 
case c; however, the Angstrom turbidity coefficient b as de- 
fined in equation (4) increases from case a to case c. This 
confusing behavior is due to the fact the AngstrOm turbidity 
coefficient is defined as the optical depth at the wavelength of 
1 mm (b = t(l mm), while the SOD curves shown in Figure 7a 
cross at a wavelength of about 0.7 mm. It is noted that atmos- 


pheric turbidity is more accurately defined by the optical 
depths at shorter wavelengths. In this investigation we choose 

P = 0.5 x (evaluated at x = 0.5 |im) (13) 

which is evaluated outside of the Chappius and other absorp- 
tion bands. The empirical value of 0.5 in equation (13) is 
chosen to make P approximately equal to Angstrom’s tur- 
bidity coefficient. As shown in Table 4, the approximation 
given in equation (13) more faithfully represents the changes 
of air pollution observed in the Mexico City atmosphere. 

Aerosol SOD is approximated using equation (7) for the 
spectral intervals free of absorption by gaseous components. 
The parameters of the approximation are given in Table 3 in 
columns 5 to 8. An examination of Table 3 demonstrates that 
equation (7) far more accurately represents the Mexico City 
observations than the AngstrOm turbidity given by equation 

(4). 

To isolate the influence of cirrus clouds, the value of SOD 
averaged from all cloudless days was subtracted from the val- 
ues measured on April 25, 1992. The result, shown in Figure 
8, is a nearly neutral spectral dependence, which is expected 
for large ice crystals (or water droplets). The absolute values 
of SOD for cirrus clouds measured on different days ranged 
approximately between 0.3 and 1.6. 

5.3. Temporal and Spatial Variations 

Temporal variations of reduced SOD are shown in Figure 9 
for the complete period of observations at three wavelengths: 
350 nm, 500 nm, and 850 nm. Variations of atmospheric 
transparency due to aerosol attenuation are very small during 
the period of observations: The curves shown in Figure 9 are 
relatively constant during this 2-month period. This demon- 
strates the fact that the atmosphere remains strongly polluted 
for extended periods of time, with no relief. Only on the 
evening of June 14, when strong winds blew some of the aero- 
sols out of the atmosphere, did the optical thicknesses de- 
crease significantly. A comparison of SOD values measured 
in the morning and evening on June 14 is shown in Figure 10. 
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Optical thickness decreased by about a factor of 3-4 during 
this 1-day period. 6 

A comparison of the SOD measured at the site of the Uni- 
versity campus with that measured in the center of the city is 
shown in Figure 11. In the center of Mexico City the SOD 
values were found to be about 50% higher in the UV than at 
the University and about 20% higher in the NIR. 

5.4. Aerosol Concentration and Aerosol Absorption 

The high values of SOD measured in Mexico City indicate 
that the concentration of aerosol particles in the atmospheric 
vertical column exceeds background levels typically found in 
other regions of the globe by an order of magnitude. Figure 
12 shows a set of theoretical calculations made assuming the 
normalized size distribution of continental aerosols suggested 
by Valley [1980], and using the optical parameters of dry at- 
mospheric aerosols in accordance with Zuev et al. [1973], In- 
tercomparison of the theoretical and observed curves suggests 
that the aerosol vertical column amount is approximately 
3x10 cmr 

Assuming that the aerosol size distribution may be repre- 
sented by the Junge form as expressed in equation (5), the se- 
ries of observations taken in Mexico City gave values of 3.25 
< v < 3.85. However, these values of v must be considered 
only representative for a smoothed distribution of aerosol par- 
ticles in the whole atmospheric column, because the ground- 
based measurements showed strong multimodal distributions. 

Figure 13 shows examples of aerosol size distributions 
measured near the ground in April 1992 using the photoelec- 
tncal counter AZ-5. As a rule, the aerosol size distributions 


in Mexico City are multimodal. Aerosol concentrations are 
measured in the range between 2x1 O' 1 and 3xl0 5 cm‘ 3 , the 
amount of optically active particles being approximately 3- 

5x10 cm'. This means that if the total vertical column con- 
centration of aerosol particles is equal to 3xl0 9 cm -2 (as es- 
timated above), then the height of a homogeneous aerosol 
layer must be about 0.6-1 km. Similar values are reported by 
Muhlia et al. [1989], 

If we assume that the part of the SOD represented by C is 
caused mainly by absorption, then a lower bound on single- 
scattering albedo EJ 0 may be estimated. Figure 14 shows 
these values of GJ 0 calculated using the average curves for 
cloudless days of the groups "a” and "c," as well as the aver- 
age curve for all the cloudless days of observations. Aerosol 
single-scatering albedos for the ensembles of aerosol particles 
observed in Mexico City have values of 0.7 - 0.9 in the UV 
portion of the spectrum, 0.6 - 0.8 in the VIS, and about 0.4 - 
0.7 in the NIR. Figure 14 also shows single-scattering al- 
bedos measured in the laboratory for room dust, coal slag, and 
soot [Lifshits et al., 1966], Note that the laboratory samples 
generally indicate an increase in single-scattering albedo with 
increasing wavelength. However, the measurements in the 
Mexico City atmosphere show the opposite wavelength de- 
pendence, because they characterize ensembles of aerosol 
particles, which scatter the light more in the shortwave range 
than in the NIR. e 

Measurements and calculations of the single-scattering al- 
bedo of aerosols over the whole spectral range for urban re- 
gions are reported in the literature. Hanel et al [1990] re- 
ports values of W Q = 0.67 - 0.76 for dry aerosols in 
rrankfurt, Germany, Ohta etal. [1990] calculated values W 0 
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Figure 14. Single-scattering albedo of the ensembles of ^erage cuwe^or the group 

^JS?e3f£ 3fSS P $ s » dust; cross, coal slag; train- 

gle, soot. 


= 0.8 - 0.9 in Sapporo, Japan, on the basis of measured 
chemical compositions. Mita and Isono [1980] estimated that 
for 30% by volume soot content in aerosols the single- 
scattering albedo may decrease to a value of 0.66. Also.Lm 
et al. [1973] found values of GJ 0 = 0.68 - 0.77 for New York 
aerosols. 

5.5. Gaseous Absorption 

The average gaseous absorption SOD may be obtained by 
subtracting the approximated (smoothed) aerosol SOD from 
the corresponding original curves. Results for the averaged 
data of the whole period of observations are shown in Figure 
15 which includes the strong absorption bands of ( 1 ) ozone 
(0 53 - 0 63 pm); (2) water vapor (0.69, 0.72, 0.82, and 0.94 
pm); and (3) oxygen (0.76 pm). Also several weak absorption 
bands are recognized, especially in the UV and VIS We have 
not identified all of them, but the maxima on the slope of the 
ozone absorption band (at wavelengths of 622 nm and 
644 nm) probably belong to NO 3 [DeMore et al, 1985], 
which is undoubtedly present in the Mexico City atmosphere. 
The curve needs more detailed and accurate investigation o 
identify all of the chemical species However, this is beyond 
the scope of the present investigation. 

6. Conclusion 

The observations in the Mexico City atmosphere confirm 
that the concentration of aerosol particles exceed by an order 


of magnitude the values found in rural regions. Analysis of 
measured values of aerosol SOD permit the concentration of 
aerosol particles to be estimated for the entire vertical column. 
On some days this value was found to be as large as 
3xl0 9 cm" 2 . 

In the UV (at X = 0.35 |im), aerosol optical thicknesses were 
found to range between 0.6 and 1 . 2 ; in the visual portion of 
the spectrum (at X = 0.5 pm) they ranged between 0.5 and 
0 . 8 ; and in the IR (at X = 0.85 pm) the optical thickness 
ranged between 0.3 and 0.5. Aerosol particles in the polluted 
atmosphere of Mexico City were strongly absorbing_ Esti- 
mates of the single-scattering albedo range between 70-90 /o 
in the UV, 60-80% in the visual portion, and to 40-70% in the 

IR 

The observations show that when the air pollution increases, 
not only do the values of SOD increase but also the spectral 
dependence of SOD changes, becoming more characteristic of 
Mie scattering. The solar radiation scattered by aerosol par- 
ticles in this polluted atmosphere is of the same order as the 
radiation scattered by air molecules. In situations of strong 
pollution the absorption by aerosols in this atmosphere can 
reach the order of the molecular absorption by atmospheric 
gases. Hence the SOD of molecular (Rayleigh) scattering, the 
SOD of molecular absorption, the SOD of aerosol scattering, 
and the SOD of aerosol absorption of light in the strongly 
polluted terrestrial atmosphere all are of comparable size. 

It was found that in the case of strong pollution the Ang- 
strom turbidity formula poorly represente the observational 
data. Applying the spectral optical depth (SOD) model of x(k 
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Figure 15, Averaged curve for gaseous absorption SOD in the Mexico City atmosphere, measured during 
April- June 1992. 



) = C + AX~ a i approximate values of 1.55 < a < 1.85 were 
obtained for polluted, cloudless days, with values of 1.25 < a 
< 1.60 on cloudless days with well-recognized haze. Almost 
neutral dependence of x versus X with valulels of 0.65 < a < 
1.05 were obtained on days with cirrus clouds. 

The Angstrom turbidity coefficient is widely used in mete- 
orological observations to characterize the degree of air pollu- 
tion. This approach appears to be adequate in atmospheres 
with weak and perhaps moderate air pollution. However, in 
cases of heavy air pollution, as found in Mexico City, it fails 
to properly represent the actual atmospheric conditions. It is 
found that one-half the value of optical depth at a wavelength 
of 1=0.5 mm is a more faithful characterization of cases with 
heavy air pollution. According to the Angstrom classification 
scheme, the Mexico City atmosphere is defined only as "very 
turbid." 

Variations of atmospheric transparency (i.e., aerosol SOD) 
were very small during the 2 months of observations; that 
means that the Mexico City atmosphere remains strongly 
polluted during extended periods of time. Also, air pollution 
in the center of the city is much higher than at the site of the 
University campus in the southern region of the city. 

Measurements of the aerosol size distributions made near 
the ground (at shelter height) showed that as a rule they were 
strongly multimodal. The observed concentrations were in the 
range between 2x*10^ and 3x10^ cm'^. This means that the 
height of the homogeneous aerosol layer must be about 0 . 6 - 1 .0 
km. 

Calculated spectral dependence of gaseous absorption shows 
not only well-known strong absorption bands of ozone, water 
vapor, and oxygen but also several weak absorption bands 


which were not identified fully. Absorption features at X = 
622 nm and X = 644 nm probably belong to NO 3 . 

Observations of SOD in Mexico City in 1992 showed with- 
out any doubt that atmosphere above the city is tremendously 
polluted. Only strong winds (e.g., June 14) significantly im- 
proved the situation; however, they are seldom found in 
Mexico City. 
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